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Abstract. Mode conversion occurs when a wave passes through a region where the sound and 
Alfven speeds are equal. At this point there is a resonance, which allows some of the incident 
wave to be converted into a different mode. We study this phenomenon in the vicinity of a two- 
dimensional, coronal null point. As a wave approaches the null it passes from low- to high-/3 plasma, 
allowing conversion to take place. We simulate this numerically by sending in a slow magnetoacous- 
tic wave from the upper boundary; as this passes through the conversion layer a fast wave can clearly 
be seen propagating ahead. Numerical simulations combined with an analytical WKB investigation 
allow us to determine and track both the incident and converted waves throughout the domain. 
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INTRODUCTION 

At the conversion region the plasma /3, the ratio of the gas pressure to the magnetic 
pressure, is approximately equal to unity. This layer generally lies low in the solar 
atmosphere, in the chromosphere. However, if a coronal null point is considered, where 
the magnetic field B — 0, there will be an area surrounding the null where /3 ~ 1 . Thus, it 
is reasonable to expect mode conversion to occur as a wave approaches the null, passing 
from low- to high-/3 plasma. 

Extensive work has been done in examining mode conversion in one-dimension. 
Zhugzhda and Dzhalilov [1] looked at the conversion of fast and slow magnetoacoustic 
waves propagating from high- to low-/3, in a gravitationally stratified, isothermal atmo- 
sphere. Transmission and conversion coefficients were found, using an exact solution. 
McDougall and Hood [2] studied conversion of slow to fast magnetoacoustic waves 
propagating downwards from low- to high-/3. Using a method developed by Cairns and 
Lashmore-Davies [3] transmission and conversion coefficients were found. This method 
has the advantage that an exact analytical solution need not be known. This was extended 
to include a non-isothermal atmosphere [4], where it was found that the same conversion 
and transmission coefficients apply. Mode conversion around a two-dimensional, mag- 
netic null point has also been investigated [5]. This looked at a fast wave propagating 
towards the null point from above, and focused on how the proximity of the /3 ~ 1 layer 
affects the competing refraction and conversion effects. 

We look at driving a slow wave pulse along the field lines from above. This travels 
from low- to high-/3 towards the null point, and conversion is observed as it crosses the 
c s = va layer. We use the WKB method to track the incoming slow wave, and as this hits 
the conversion region we also track the fast wave which propagates ahead. 
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FIGURE 1. Equilibrium magnetic field with a magnetic null point situated at the origin, denoted by a 
cross. The circle shows where the sound and Alfven speeds are equal. 



MODEL 



We use the ideal MHD equations with gravity neglected. The equilibrium magnetic field 
is given by Bo = 5o(x,0, — z) /L, and is shown in Figure 1. The squared sound and 
Alfven speeds are given by c 2 = YPo/Po an d v 2 A = Bq (x 2 +z 2 ) / (/ipo-L 2 ) respectively, 
where 7 is the ratio of specific heats, po and po are the constant equilibrium pressure and 
density, /1 is the magnetic permeability, and L is a typical coronal length scale. 

The ideal MHD equations may be linearised and combined to give a pair of linear 
wave equations: 
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We define the velocity parallel and perpendicular to the magnetic field as 

t2 respectively. We may then 
and vii = sin cot, 
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V|| = (xv x — zv z ) /y/x 2 +z 2 and = (zv x +xv z ) 
drive a slow wave pulse exactly along the field lines by setting vj 
where 00 is the driving frequency. As the pulse is curved compared to our horizontal 
boundary, we calculate when each point hits the boundary, so the pulse is released along 
the boundary at the correct time. 



WKB METHOD AND SIMULATIONS 

Starting with the wave equations (1) and (2), we write v x = aexp(?0 (x, z, t)) and 
v z = bexp (i<j> (x, z, t)) under the assumption that <j> ^> 1. We then obtain the quadratic 
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which we have set equal to # ', where p = d(j)/dx and q = d§ jdz- The roots of this 
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FIGURE 2. Contour plots of the parallel velocity from the numerical simulation. The circle denotes 
where c s = Va, and the cross is the location of the null point. The lines show the front, middle and back of 
the slow wave and fast wave as predicted by the WKB method. 



quadratic are 
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where the plus sign gives the fast wave solution, and the minus sign the slow wave. 

Using Charpit's relations we may then obtain a system of ODEs governing the fast 
and slow magnetoacoustic waves: 
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We solve these using a fourth order Runge-Kutta scheme, allowing the position of the 
wave pulses to be tracked throughout the domain, as shown in Figure 2. 

The slow wave is tracked as it enters the domain and approaches the conversion layer. 
Once the front of the pulse reaches this point, we begin to follow the fast wave pulse 
that is created. This is only tracked for points which pass through the conversion layer, 
hence a much smaller portion of the wave front is followed. 

CONCLUSIONS 

We have studied a two-dimensional mode conversion layer situated around a magnetic 
null point. A slow wave pulse is driven along the field lines on the upper boundary, 
and as it hits the conversion region, where c s = v A , some of its energy is transferred to 
the fast mode, which we see propagating out ahead of the slow mode pulse. The WKB 
method is then used to predict the position of both wave fronts as they travel through 
the domain. These predicted positions are in excellent agreement with the simulations, 
showing the slow wave stretching out and slowing as it approaches the null, while the 
fast wave propagates out in front. 

In future we plan to extend the Cairns and Lashmore-Davies [3] method to two 
dimensions, allowing the quantity of conversion and transmission to be calculated, as 
has been done for an isothermal [2], and non-isothermal [4], one-dimensional model. 
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